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ABSTRACT 

We present simultaneous high-resolution optical spectroscopy and X-ray data of the 
X-ray binary system GRMus (XB 1254-690), obtained over a full range of orbital 
phases. The X-ray observations are used to re-establish the orbital ephemeris for this 
source. The optical data includes the first spectroscopic detection of the donor star 
in this system, through the use of the Doppler Tomography technique on the Bowen 
fluorescence blend (~4630-4650 A). In combination with an estimate for the orbital 
parameters of the compact object using the wings of the Hen A4686 emission line, 
dynamical mass constraints of 1.20 ^ Mx/Mq ^ 2.64 for the neutron star and 0.45 
^ M2/MQ ^ 0.85 for the companion are derived. 

Key words: binaries: close - stars: binaries: individual: GRMus, XB 1254-690 - 
binaries: spectroscopic - X-rays: binaries. 



1 INTRODUCTION 

XB 1254-690 is a persistently bright, low-mass X-ray binary 
(LMXB). It was identified with a faint blue star (GRMus, 
V — 19.1), whic h exhibited the Bow en blend of NIII and 
C III in emission ( Griffiths et al ]| 19781). The object produced 

Hourvoisier et al.l 



Type 1 X-ray bursts l|Mason et alj|l98(il [O 



1 19861 ). indicating the presence of a neutron star as the com- 
pact object. Dips of up to 95% of the 1-10 keV flux, with 
a recurrence period of Px = 3.88 ± 0.15 hr and lasting 
~ 0.8 hr per cyc l e, we re discovered in EXOSAT data by 
ICourvoisier et all fl986). The dips are caused by obscura- 
tion of the central source by a bulge on the outer edge of 
the accretion disc, implying a moderately high inclination. 

From F - band observations of the optical counterpart, 
iMotch et all j 19871 ) determined an optical ephemeris where 
minimum light occurs at T mi „ = JD 2,445,735.693 (4) + 
0.163890 (9) days, equivalent to a period of P op t = 3.9334 
± 0.0002 hr. This broad optical modulation is probably due 
to the changing visibility of the heated face of the secondary 
star in the system, with an additional contribution from the 
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X-ray heated bulge, although other explanations such as an 
asymmetrical disc or occultation of the disc by the com- 
panion cannot be ruled out. Within the errors, the optical 
and X-ray period measurements are consistent, with later 
X-ray observations indicating that the optical modulation 
provides a more accurate measurement of the orbital period 
jSmale et alj|2002l ; iLevine fc Corbetll2006l ). A short section 
of simultaneous X-ray and optical coverage showed that the 
optical minimu m occurs ~ 0.16 in phase after the centre of 
the X-ray dips (|Motch et al.lll987f ). 

Simple geometrical modelling is possible by considering 
the X-ray illuminated area of the companion star, disc ra- 
dius, disc opening angle, and inclination at given distances 
which are compatible with the observed optical flux and the 
amplitude of the light curve. This leads to constraints on 
the sourc e inclination of 65 °-73°, and on the distance of 
8-15 kpc (|Motch et al.ll 19871 ). 

The first mass constraints for this system were based 
upon a unique relation between the orbital period and the 
mass of the companion, assuming that the companion is a 
low mass, zero - age main sequence star fill ing its Roche lobe 
l|Warnerll 19761 ). Using this approximation. ICourvoisier et al.l 
(1986) estimate the mass of the companion to be in the 
region of ~ 0.45Mq. 
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Table 1. Log of observations of GRMus / XB 1254-690 



HJD 


Observatory 


Total Int. Time 


2453151.52 


VLT 


3.75 h 


2453151.51 


RXTE 


4.75 h 


2453152.52 


VLT 


4 h 


2453152.69 


RXTE 


5.25 h 


2453153.53 


VLT 


2 h 



The first kinem atic mass constraints were proposed by 
ICowlev et al] |l988), who obtained intermediate resolution 
(~ 3A) spectra using the 4m telescope at the Cerro Tololo 
Inter-American Observatory. Assuming that velocities de- 
rived from the Hell A4686 emission line centroid represent 
the motion of the compact star and not, for example, the 
streaming of gas between the stars or out of the system, then 
approximate stellar masses can be derived using the mass 
function. The result is consistent with a low-mass main- 
sequence star (M ~ 0.5-0.8 Mq) and a canonical 1.4 Mq 
neutron star. 

In the remainder of this paper we will attempt to con- 
firm and better constrain these mass estimates using high- 
resolution blue optical spectroscopy. We will measure the 
wings of the He II A4686 emission line to gain an estimate for 
the compact object veloci ty (Section 3.3), and use t he Bowen 
Fluorescence technique (Ste eghs fc Casaresl I2002T ) to mea- 
sure the donor star velocity (Section 3.4). We also obtained 
simultaneous X-ray observations (see Section 3.1) in order 
to re-establish the ephemeris of iMotch et al] (|l987f) . This 
cannot be extended to the observations described herein be- 
cause the propagated uncertainty in phase is now compara- 
ble to the orbital period itself. 



2 OBSERVATIONS 

We have obtained simultaneous optical and X-ray data for 
GR Mus / XB 1254-690. The journal of observations is given 
in Tabled] 

2.1 Optical data 

Phase- resolved spectroscopy of GRMus was obtained from 
26 to 28 May 2004 using the FORS 2 spectrograph mounted 
on the VLT/UT4 at ESO's Paranal observatory. During each 
night we observed GR Mus for ~ one full orbit, resulting in a 
total of 38 spectra with an integration time of 900s each. We 
used the 1400V volume-phased holographic grism, and a slit 
width of 0.7", giving a wavelength coverage of AA4513-5814 
with a resolution of 70 kms" 1 (FWHM). The seeing during 
these observations varied between 0.5 and 2.1 arcsec. The 
slit was orientated at a position angle of 88.15° to include a 
comparison star in order to calibrate slit losses. We observed 
the flux standard Feige 110 with the same instrumental set- 
up in order to correct for the instrumental response. He, 
Ne, Hg and Cd arc lamp exposures were taken during day- 
time for the wavelength calibration scale. We de-biased and 
flat-fielded all of the images and used optimal extraction 
techniques to m aximise the signal-to-noise ratio of the ex- 
tracted spectra l|Horneil 19861 ) . The pixel-to- wavelength scale 
was derived through polynomial fits to a large number of 
identified reference lines resulting in a dispersion of 0.64 A 



pix . Velocity drifts due to instrumental flexure (found to 
be always <5 kms -1 ) were corrected by cross-correlation 
with the sky spectra. 



2.2 X-ray data 

Simultaneous X-ray data from XB 1254-690 were obtained 
using the Proportional Counter Array (PCA) on board 
the Rossi X-Ray Timing Explorer (RXTE). The Propor- 
tional Counter Arr ay (PCA; for a detailed description, see 
IJahoda et all (|l996l )) onboard the R XTE satellite consists of 
an array of 5 co-aligned Proportional Counter Units (PCUs) 
that are sensitive to photons of energy 2 to 60 keV with a 
total collecting area of 6500 cm 2 . We obtained 36 ks of total 
coverage during our optical observations. For our analysis 
we only used the data from the Standard 2 configuration, 
with a time resolution of 16 seconds. 



3 ANALYSIS 

3.1 Determination of orbital phase 

The ephemeris of IMotch et al] (|l987l ) is now out of date be- 
cause the propagated uncertainty in phase is over a factor of 
2 greater than the orbital period. We therefore need to make 
use of our contemporaneous RXTE observations to set a new 
zero point u s ing th e X-ray dips (note that we still use the 
IMotch et al] (| 19871 ) orbital period) . This is not quite as sim- 
ple as it would first appear since there are large variations in 
dip depth and structure that are apparent from observation 
to ob servation, as well as cycle to cycle l|Courvoisier et all 
1986). Additionally, a note of caution must be sounded since 
in a number of earlier RXTE and BeppoS AX observations 
l|Smale fc Wachterl Il999l ; llaria et al.ll200ll ). X-ray dipping 
was not observed at all, with an accompanying decrease in 
the mean optical variability of the source. With the com- 
bination of these effects, it was suggested that the vertical 
structure on the accretion disc edge, associated with the im- 
pact point of the accretion stream, had decreased in angu- 
lar size at the time of these earlier observations. The cause 
of this intermittent reduction in opening angle is not un- 
derstood. Our observations include just one deep dipping 
episode (Fig. [TJ, co- incident with our first night of optical 
observations. 

XB 1254-690 typically exhibits 'shoulders' of dipping, 
i.e. regions both before and after the deep dipping in which 
there is a small decrease of intensity. Each shoulder of the 
dip lasts ~650-925 seconds, during which the total count 
rate in the 2-25 keV band drops by ~15% with no strong 
changes in the hardness ratio. The shoulder is followed by 
an inter val of strongly vari able deep dipping lasting ~1100 
seconds (|Smale et al.ll20o"2l ). with associated strong changes 
in the hardness ratio (Fig. [2j . The hardening of the spectra 
in the dips is most likely to be due to absorption by cold 
matter. 

Our observation of a dipping episode opens with the 
system already in the first shoulder phase, the onset of deep 
dipping occurring approximately 200-250 seconds into the 
observation (before a brief return of the source intensity to 
the non-dip level), and then continuing to capture the full 
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Figure 1. 2-20 keV X-ray lightcurve showing the deep dip used 
to calibrate a new ephemeris for the system. The dipping episode 
appears to begin shortly before JD 2453 151.64, and occurs for the 
typical dip duration of ~1100 seconds jSmale et alj|2002h . 
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Figure 2. 2-4 keV (top) and 4—9 keV (centre) X-ray lightcurves 
displaying the deep dipping episode. The bottom panel displays 
hardness as a function of intensity. There is an anti-correlation 
between intensity and hardness ratio during the deep dipping 
episode, which begins shortly after 200 seconds on the above scale. 
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Orbital Phase 

Figure 3. Continuum lightcurve of GRMus, normalised to the 
flux level of the first observation in our dataset. Note that the 
errors on each measurement are smaller than the plotted points. 

the dip centre by eye, taking the mid-point between the on- 
set of and exit from deep dipping (which can be determined 
with the help of the associated hardness ratio changes, see 
Fig. H>. This occu rs at T dip = JD 2453151.647 (3). From 
iMotch et al l (|l987t ) we know that the dip centre has an av- 
erage orbital phase of 0.84 (with phase occurring at optical 
minimum). We therefore estimate a new value for To of JD 
2453151.509 (3). 

There are problems with the reliability of an ephemeris 
generated from just a single dipping episode, not least be- 
cause we do not have full coverage either side of the dip. It is 
therefore hard to estimate the dip centre with absolute con- 
fidence. However, we are encouraged by the fact that our dip 
lies within 0.05 phase of the 1999 dips (the propagated phase 
un certainty from these d ips is of the order ~0.6) reported 
bv lin't Zand et all (|2003l ) . More importantly, the continuum 
lightcurve from the resulting ephemeris possesse s the same 
phasing relation reported by IMotch et all (|l987l ) (Fig. [HjEl- 
In light of this evidence we are confident that we have iden- 
tified the dip centre correctly and we will use this zero-point 
to calculate orbital phases throughout the remainder of this 
paper. 



egress from deep dipping, into the shoulder and then re- 
turning to the continuum level. Note that the return of the 
source intensity to the non-dip level during this episode is 
not unusual, and the light curve can in f act be interprete d 
as superpositions of small micro-dips (e.g. lUno et alj|l997l ). 

Unfortunately, the lack of information regarding the full 
ingress to deep dipping through the preceding shoulder in- 
troduces a degree of uncertainty into our estimate of the dip 
centre. For example we could potentially be observing only 
a small section of an extended deep dipping episode (e.g. 
lUno et al1ll997l . report dips lasting for ~40 minutes). If this 
is the case, the absence of the beginning of the dip will cause 
the estimate for dip centre to be incorrect. Nevertheless, we 
do have cause for optimism. The deep dipping which we ob- 
serve is approximately 1100 seconds in length, equivalent 
to the typical period of a dipping episode in XB 1254-690 
(see above). This is a strong indicator that we have full, or 
nearly-full, coverage of the deep dipping phase. We estimate 



3.2 Optical emission spectrum 

The continuum subtracted average of the 38 individual op- 
tical spectra, corresponding to an exposure time of nearly 
10 hours, is given in Fig. [4] The overall spectrum resembles 
that of other LMXBs, with high excitation lines superposed 
on a blue continuum. This is typical of strongly X-ray h eated 
accretion discs (e.g. Ivan Paradiis fc McClintocklll995l ). 

The blue end of the optical spectrum in this system 
is dominated by the Bowen emission blend of C ill and 
Nni. Hell A4686 is also seen strongly in emission, in ad- 
dition to weaker H/3 and Hell A5411. The absorption com- 
ponent longward of H/3 is observed in a few other LMXBs 

1 Note that the large intrinsic scatter in the optical continuum 
lightcurve makes the generation of a new ephemeris solely from a 
relatively limited dataset such as ours rather difficult. 
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Figure 4. Average normalised spectrum of GRMus. The most 
prominent emission lines are indicated, including the Bowen blend 
of Cm and Nm, Hen and H/3. The emission line near 5577A is a 
night sky feature. 



such as X 1822-37 1 dCasares et alj l2003h. XTE J2123 -058 
ilHvnes et al.ll200ll). GX339 -4 l|Buxton fc Vennesll2003h and 
MMSer (|Hvnes et al.|[20oj ). 

The narrow emission feature at 5577A is [Ol], a strong 
night-sky line useful for confirming the veracity of our wave- 
length calibration. The broad features between 4900-5300A 
remain unidentified. The absence of these features in our 
comparison star lead us to believe that they are real. We also 
see some of these 'bumps' in other LMXBs (e.g. GX9+9, 
X 1822-371, V801Ara, V926Sco), and perhaps they are 
blends of a number of different weak emission lines. Part 
of this could be He I A4922 and A5016, and also Fell A4924, 
A5018, A5169 and A5297. There are also several O n and 
O in lines between 4891-5006A which may explain the large 
bump at ~4930A. Note also that the Bowen blend (and to a 
lesser extent the He II A4686 emission line) appears to be sit- 
ting upon a broad emission 'h ump'. This has been noted in 



other LMXBs (e.g. V801 Ara; lAugusteiin et al.lll998T). and 
attrib uted to a blend of Fe II emission lines ( Schachter et al.l 
1 19891 ). 

In Fig. [5] we display the Bowen blend and He n A4686 
in the form of a trailed spectrogram. All of the individual 
spectra were first normalized to the continuum using a third- 
order polynomial fit, and then phase- folded into 15 bins. The 
spectra were then smoothed using a Gaussian filter with a 
FWHM of two pixels. We see a clear S-wave structure mod- 
ulated on the orbital period in the He II A4686 emission. The 
blue-to-red crossing phase is ~0.6 and the semi- amplitude 
~250 kms -1 . This behaviour indicates a close association 
with the accretion disc or compact object. There is also an 
additional, weaker, component present. Again, this is mod- 
ulated on the orbital period, approximately a quarter cycle 
out of phase with the more prominent component. There 
also appear to be faint traces of two anti-phased sharp com- 
ponents in the Bowen blend though the broad underlying 
component moves in phase with the accretion disc S-wave 
(maximum blue-shift at phase 0.25). The tomograms dis- 
cussed in Section 3.4 will help to characterise these spectral 
trails more thoroughly. 

The H/3 line is shown in the right-hand panel in Fig. [5] 




Velocity (km s ) 

Figure 6. Velocity profile of the Hen A4686 emission line dis- 
playing a complex core structure. The extreme blue wing is con- 
taminated by the Bowen emission blend. 



along with the broad emission feature at ~4930A. The ab- 
sorption component longward of H/3 is present at all phases, 
and in common with the emission component displays an 
S-wave approximately in phase with the expected motion of 
the accretion disc. The broad feature at ~4930A also dis- 
plays a S-wave. The Hell A5411 emission line is too weak in 
the individual spectra to produce a trailed spectrogram. 



3.3 Motion of the compact object 

Early attempts to determine any kine matic constra i nts on 
this system met with limited success. iMotch et all (|l987f ) 
noted that the relative weakness of the emission lines in 
their spectra did not allow for the study of their variation 
with orbital phase. However, they did exclude the possibil- 
ity of orbital motion of He II A4686 (which requires a source 
of highly ionising radiation, and is typically seen in accre- 
tion discs) in exces s of 300 kms" 1 . Later observations by 
ICowlev et alj l|l98ct ) did result in a measurement of the same 
Hell line's core velocity, which possessed a semi-amplitude 
of 114 ± 13 kms -1 . However, this is not a particularly re- 
liable measurement of the compact object velocity. The line 
core can often be heavily contaminated by other emission 
sites (see Fig. [SJ), for example from the accretion stream or 
'hot-spots' on the accretion disc itself. 

Instead, the compact object velocity can be better es- 
timated by examining the wings of the emission line. These 
wings find their origin in the high- velocity gas of the inner 
accretion disc, very close to the compact object itself. Pre- 
sumably this gas will closely replicate the motion of the com- 
pact object, and it is far less likely that alternative emission 
sites will provide a strong contaminating effect. However, 
this method is not immune to contamination and we must 
be cautious in accepting the results of this analysis. 

We have followed the double-Gaussian technique 
l|Schneider fc Yound ll980T ) to estimate both the velocity 
semi-amplitude K\ and the systemic velocity 7 of the com- 
pact object from Hell A4686. We employed a double Gaus- 
sian bandpass with FWHM of 100 kms -1 and Gaussian 
separations of 300-1700 kms -1 in steps of 100 kms -1 . Ra- 
dial velocity curves of line sections with separations of 700- 
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Figure 5. Trailed spectra of the Bowen complex and Hell A4686 (left) and the H/3 region (right). The data are plotted twice for clarity. 
Note that the wavelength scales differ. There are clear multiple S-waves evident in the Hell emission line, ~0.25 out of phase with each 
other, and also in the Bowen complex. The absorption component redwards of H/3 is present at all phases. 



1400 kms -1 yield consistent results with a median blue- 
to-red crossing phase of 0.58 ± 0.03, median velocity semi- 
amplitude K\ = 130 ± 16 kms -1 and median systemic ve- 
locity 7 = 173 ± 12 kms -1 (Fig. [7J. The errors are conser- 
vatively taken from the maximum error obtained from sine 
wave fits to the results from each individual line section. 
Separations below 700 kms -1 suffer badly from contamina- 
tion by the line core, whilst separations over 1400 kms -1 be- 
come corrupted by continuum noise and the Bowen emission 
blend. The small difference in phasing (0.08 orbits; slightly 
larger than the ephemeris uncertainty of ~0.02 orbits) in- 
dicates that even the line wings may be contaminated to 
some extent by inhomogeneities in the disc. This contam- 
ination will also affect the measurements of K\ and 7 to 
some extent. Nevertheless, this remains the best method for 
estimating K\ from our own and other currently available 
data. 



Unfortunately, the Hell A5411 and H/3 emission lines, 
which are clearly seen in our averaged spectrum (Fig. [4}, 
are too weak in the individual spectra for a similar analysis. 
The fitting procedures are simply dominated by noise. 



3.4 Doppler tomography 



We can use Doppler tomography l|Marsh fc Hornelll988l ) to 
measure the radial velocity amplitudes of any weak, narrow, 
emission lines that may be present. By mapping the observed 
data onto a velocity coordinate frame, Doppler tomography 
makes use of all data at once and can thus be used to search 
for features that are too weak to be distinguished in each 
individual spectrum. One effectively resolves the distribu- 
tion of line emission in the co-rotating frame of the binary 
system, providing an excellent tool for identifying the ori- 
gin and kinematics of the various emission components. It 
also allows us to map the emission line distribution from the 
broader He 11 and H/3 lines. 

Secondary star emission is readily identifiable in 
Doppler tomograms since the solid body rotation of the 
Roche lobe is mapped to a corresponding Roche lobe area 
along the positive Vy-axis. For reference, in Figs. 151 4111 we 
overplot the area of the Roche lobe and the centres of 
mass of the donor star (uppermost cross), full system (map 
centre) and the compact object (lower cross) assuming a 
ma ss ratio of ~0-4 (eq uivalent to the mass ratio proposed 
by ICowlev et a.1.1 [T988f ) and compact object velocity of 130 
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Figure 7. Diagnostic diagram for the Hell A4686 emission line 
wings, displaying the systemic velocity 7 (top), phasing (centre) 
and radial velocity K\ (bottom). 




kms' 1 . Additionally, we have overplotted the ballistic ac- 
cretion stream for such a system. The circle, centred upon 
the compact object and with radius 600 kms -1 , indicates 
the common emission regions from a generic LMXB accre- 
tion disc. We employed a maximum entropy implementation 
of Doppler tomography whereby the data are fitted under 
the added constraint of maxim izing the entropy (i.e. image 
smoothness) of the tomogram (|Marsh fc Hornelll988l ). This 
reduces the presence of noise artifacts in the recovered to- 
mograms and allows for a simultaneous fit to a number of 
heavily blended lines. 

In Fig. |S]we display a tomogram of the He n A4686 emis- 
sion line. This displays a clear enhancement in the lower left 
quadrant and an additional component close to the accretion 
stream. We do not observe the classic ring-like distribution 
of the accretion disc, presumably due to the strength of the 
line core features vastly overwhelming the high velocity inner 
disc emission. For this reason we cannot be fully confident 
in the results of the double Gaussian technique (see Section 
3.3), though this remains at present the best method for 
estimating K%. 

The enhancement in the (lower) left quadrant is also 
seen in the tomograms of H/9 (Fig. [9} and Hen A5411 (Fig. 
UOp . Whilst this phenomenon is a common feat ure in the re- 
markable SW Sex-type systems l|Hellier|[2OO0l ). it has also 
bee n noted in some other LM XBs, notably XTEJ2 123- 
058 |Hvnes et al.ll200lf), GX9+ 9 (jCornelisse et al.l2006l ) and 
J 0422+32 l|Casares et al.ll 19951 ). However, in this case (and 
in the case of X 1636-536 and X 1735-444, ICasares et all 
2006) the emission is not quite so localised in the lower- left 
quadrant, but is present in a broad band on the left-hand 
side of the maps. Rather than being an analogue to the 
highly magnetic SW Sex-type systems, it is likely that this 
emission is produced in an extended disc bulge. In most 
LMXBs, the Hen emission is predominantly produced in a 
hotspot where the accretion stream impacts the disc (thus 
producing a spot in the uppe r left quadrant of a Doppler 
map, e.g. IPearson et aJj |2006) or is more evenly distributed 
around the bulk of the accretion disc. If the emitting gas was 
symmetrically distributed around the primary, the Doppler 
image would reveal a circular distribution centred on V x =0, 







Velocity (km/s) 



500 



Figure 8. Doppler tomogram reconstructed from the Hell A4686 
emission line using maximum entropy optimisation, 7 = 185 
kms -1 . The area of the Roche lobe, the ballistic accretion stream 
and the centres of mass of the donor star (uppermost cross) , sys- 
tem (map centre) and the compact object (lower cross) are plot- 
ted assuming a mass ratio of ~0.4 and compact object velocity of 
130 kms -1 . The circle is centred upon the compact object and 
possesses a radius of 600 kms -1 . 
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-K\ . This would have been a useful way of confirming 



the K\ estimate derived from the wings of the Hell A4686 
emission line. 

Like Hell A4686, the Hell A5411 tomogram also in- 
dicates the presence of some diffuse emission at velocities 
which partially correspond to the Roche lobe of the donor 
star. However, as this emission is diffuse in the Doppler 
map, extending over a broad range of velocities, it is dif- 
ficult to directly associate this extended emission with the 
donor star. Fluorescent emission from the donor in this case 
should manifest itself as a tightly concentrated spot situated 
on the positive V y axis. 

In Fig. [11] we present a tomogram of the N ill A4640 
emission line. The maximum entropy method of producing 
Doppler tomograms allows for a simultaneous fit to a num- 
ber of different lines to produce one map. However, in this 
case the N ill A4640 line strongly dominates over any oth- 
ers in the Bowen blend. We see two regions of enhanced 
emission. The brighter of the two spots occurs at a com- 
bination of phasing and velocity where we would expect 
the donor star to be fluorescing. The small horizontal off- 
set of Vi~35±30 kms -1 is consistent with the errors in our 
ephemeris (± 0.02 in phase). 

The second spot occurs at a much lower velocity, and is 
1-2 a weaker (relative to the local noise level) than the first. 
It is difficult to provide a truly satisfactory explanation for 
its origin. Perhaps it traces some form of low velocity out- 
flow (e.g. disc wind) , though why it should be preferentially 
observed at an apparent phase of ~0.25 remains a mystery 
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Figure 9. Doppler tomogram reconstructed from the H/3 emis- 
sion line using maximum entropy optimisation, 7 = 185 kms . 
The symbols are as for Fig. [8] 
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Figure 10. Doppler tomogram reconstructed from the Hen 
A5411 emission line using maximum entropy optimisation, 7 = 
185 kms . The symbols are as for Fig. [8] 
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Figure 11. Doppler tomogram reconstructed from the N III A4640 
emission line using maximum entropy optimisation, 7 = 185 
kms -1 . The symbols are as for Fig. [8] 



(orbital phase can be considered to increase in a clockwise- 
fashion around a Doppler map). 

If we presume that the bright spot in Fig. [TT] is caused 
by the irradiated inner face of the donor star, this is an excel- 
lent tracer for the kinematics of the system. The precise lo- 
cation of the spot can be measured using a two-dimensional 
Gaussian fit (though bearing in mind that the actual shape 
of the spot is not expected to be Gaussian) , providing a ve- 
locity of Kcm = 245 ± 30 kms -1 . Using this value we can 
shift all of our spectra into the rest-frame of the (irradiated 
face of the) donor. A sharp N ill A4640 component is clearly 
resolved (Fig. I12|) . However, since this is a measure of the 
motion of the irradiated inner face of the donor star, it only 
provides us wi th a lower-limit on the velocity of the centre of 
mass (see e.g. ISteeghs fc Casaresll2002j ; |M ufioz-Darias et al.l 
l2005h . 

One final piece of information remains to be extracted 
from the Doppler tomograms; assigning the correct systemic 
velocity should produce the most sharply defined features in 
a Doppler map. We used the Hell and Nm Doppler maps 
to verify the systemic velocity. The \ 2 value of the map 
was calculated for a range of 7, and the best fit in terms 
of minimal \ 2 was achieved for 7 = 192 ± 2 km s" 1 for 
the Nm map and 7 = 185 ± 2 km s -1 for Hell A4686. 
These are both consistent with, if slightly higher than, the 
results derived from the Hell A4686 wings in Section 3.3. 
The complex, asymmetric core profile of the Hell emission 
line could be contaminating one wing more than the other, 
leading to this small offset. We heretofore adopt the value for 
7 derived from the Hell A4686 Doppler map. This should be 
more reliable than the N ill map since only a single, isolated 
emission line is involved in the fitting process, rather than 
the more complicated blend of Bowen lines. 
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Figure 12. Average spectrum created by combining all of the 
individual unshifted spectra from our observations (top) , showing 
the Bowen emission blend and He II A4686 region. Below this is 
shown the average spectrum when all of the individual spectra 
are shifted into the rest-frame of the donor star. A sharp Nni 
A4640 component is clearly resolved. 

4 DISCUSSION 

4.1 Systemic velocity 

Whilst the systemic velocities inferred from both the wings 
of the Hen A4686 line and the Doppler tomograms are 
self-consistent (~ 185 kms - ), they are considera bly higher 
than the earlier estimate bv lCowlev et al.l l| 19881 ) (98 ± 10 
kms -1 ). This measurement was obtained from rather low 
quality data, with poor spectral and temporal resolution. 
Additionally, they based their measurements on the motion 
of the He II A4686 line core. There are a number of different 
emitting regions which contribute to this line, leading to a 
complex core structure (see e.g. Figs. . A basic mea- 
surement of the line centre under the assumption that this 
will accurately trace the motion of the compact object will 
therefore be subject to large systematic errors. 

In light of this, we believe that our estimate for the 
systemic velocity will be a better reflection of the true value. 
Although large, our measurements do seem reasonable in the 
context of the high galactic latitude (b = -6.42) of GRMus; 
the neutron star is likely to have received a substantial 'kick' 
out of the Galactic Plane at birth. A number of statistical 
studies on pulsar velocities have been carried out. These 
studies give a mean birth velocity of 100-500 k m s~ , with 
some possessing velocities of ^1000 km s _1 l|Wang et al.l 
l2006h . 

4.2 Interpretation of the Bowen blend Doppler 
map 

The N ill Doppler map (Fig. Hip appears to display enhanced 
emission in the region of and at the phasing where we would 
expect to see the companion star Roche lobe, in addition to 
an abundance of noisy structure in the outer regions of the 
map. A useful way to confirm the validity of the Doppler 
map is to observe the behaviour of sharp components in our 
trailed spectrogram (Fig. [5| . This is not trivial due to the 



Predicted Observed 




Figure 13. The right hand panel displays the observed trailed 
spectra of the Hell A4686 emission line and the left hand 
panel shows the ideal trailed spectrogram reconstructed from the 
Doppler tomogram shown in Fig. [8] Orbital phase is plotted along 
the x-axis. 




Figure 14. The right hand panel displays the observed trailed 
spectra of the N III A4640 emission line, part of the Bowen emis- 
sion complex. Next to this is plotted the ideal trailed spectro- 
gram reconstructed from the Doppler tomogram in Fig. 1111 Or- 
bital phase is plotted along the x-axis. 



faintness of GRMus compared to a system such as ScoX-1, 
which clearly feature s narrow Bowen lines mo ving in phase 
with the donor star (jSteeghs fc Casaresll2002l ). 

We can reconstruct an ideal trailed spectrogram from 
the Doppler tomograms to test for systematic residuals in 
the maps. As a test, we first apply this technique to the Hell 
A4686 emission line. The two distinct components we see in 
the Doppler map (Fig. [SJ are reproduced clearly in the pre- 
dicted trailed spectrogram (left hand panel, (Fig. I13[l . and 
can be easily identified in the observed data. For the case 
of the N ill A4640 emission line (Fig. I14J1 we again see two 
components in the reconstructed spectrogram, though not 
quite so distinctly as for the case of Hell. Comparing this 
idealised image to the real data it is indeed possible to detect 
faint traces of these two narrow, approximately anti-phased 
lines which produce the two distinct spots in Fig. 1111 The 
second spot is therefore unlikely to be a systematic residual 
produced when creating the Doppler map, and its unusual 
phasing/positioning could perhaps therefore be attributed 
to a region which violates one of the axioms of Doppler To- 
mography, for example that all motion should be parallel to 
the orbital plane. Note that a violation of these axioms does 
not mean that the map is erroneous, but it does complicate 
the successful interpretation of the map. 
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Table 2. Derived system parameters for GRMus / XB 1254-690 



Parameter 




Parameter 






T dlp (HJD) 


2453151.647 (3) 


-y (km s — 1 ) 


185 


± 2 


T (HJD) 


2453151.509 (3) 


9 


0.32- 


-0.43 


K% (kins -1 ) 


130 ± 16 


M X (M Q ) 


1.20- 


-2.64 


K om (kins -1 ) 


245 ± 30 


M 2 (M ) 


0.45- 


-0.85 



Note: Inclination constrained to 65°-73°: Motch et a], f 1987 1 



4.3 Refined mass limits 

Ide Jong et all (ll996T ) estimated an opening (semi-) angle for 
the accretion disc in GRMus of ~ 12°, using a simple ge- 
ometric model for the reprocessing of X-rays in LMXBs. It 
is important to note that this value is hostage to the im- 
plicit assumptions underlying these models, for example the 
geometric description of the disc and of the X-ray emitter. 
This will entail a systematic uncertainty that likely exceeds 
the internal model accuracy of le ss than a degree. N ever- 
theless, this value agrees well with iMotch et all(|l987f ) who 
estimated an opening angle of ~ 9° — 13° in order to explain 
the observed amplitude of the orbital photometric modula- 
tion. 

We can use this information to derive a mass estimate 
for the individual comp onents by using the ' K-corr ection' 
algorithm developed by iMunoz- Darias et al] (|2005l ). This 
measures the deviation between the reprocessed light centre 
of the fluorescence emission lines and the centre of mass of a 
Roche lobe filling star in a persistent LMXB, including the 
screening effects by a flared acc retion disk. Using the fourth 
order polynomial fits given by iMufioz-Darias et al.l (2005) 
for disc opening angles of 8°, 10°, 12° and 14° gives a total 
mass range for the neutron star of 1.20 ^ Mx/Mq ^ 2.64. 
Since the neutron star mass and the mass ratio are related 
by the mass functiorQ, we obtain 0.32 ^ q ^ 0.43, and limits 
on the companion mass of 0.45 ^ M2/MQ ^ 0.85. 

This correction is dependent upon the opening angle of 
the accretion disc, monotonically decreasing with opening 
angle such that the largest opening angle gives the lowest 
neutron star mass. Since the larges t opening an g le we have 
used exceeds the est i mates of both IMotch et alj (|l9871 ) and 
IMufioz-Darias et al.l (|2005l ). this gives a very conservative 
lower limit to the neutron star mass, and the true value is un- 
likely to be so low. Equally, the upper limit is derived using 
an op ening angle lower than that suggested bv lMotch et al.l 
(| 19871 ) and as such the neutron star is unlikely to be quite 
so massive. 

Limiting ourselv es to a disc opening angle of exactly 12° 
l|de Jong et alj|l996h . the mass constraints upon the neutron 
star tighten up to 1.35 «S M x /M & < 2.32, with 0.33 < q < 
0.41 and 0.49 < M 2 /Mq < 0.79. These limits are not a vast 
improvement over the ones derived above using a range of 
opening angles from 8-14°. This indicates that the major 
source of uncertainty in this case is still incurred from the 
measurements of -Kcm and particularly of K\ (see Section 
3.3), in addition to the system inclination angle. The choice 
of disc opening angle (within a reasonable range) has only 
a small effect in comparison. Nevertheless, we feel it is pru- 
dent not to restrict our final mass estimate by the assump- 

2 f(My) = K l P °^ = M X Bin 3 i 



tion of a disc opening angle of exactly 12°, preferring the 
more conservative range of 8-14°. The derived parameters 
for GRMus are listed in Tabled 

4.4 Roche lobe size 

We know that the companion must be filling its Roche lobe 
in order to permit persistent accretion. There is a well known 
and useful relationship bet ween the fractio nal Roche lobe 
radius and the mass-ratio |Paczvnskilll97ll ) which can be 
combined with Kepler's 3rd law to provide an estimate of 
the size of the Roche lobe as a function of period. In the case 

of this system (with a 3.9 hour period), the relationship may 

1/3 

be expressed as: Rl = 0.58Af 2 . 

We can thus derive the range of potential Roche lobe 
radii consistent with the mass range of the companion star 
that we have derived kinematically. This turns out to be 
0.44 < Rl/Rq < 0.55. Whilst at the lower end of our 
mass constraints the donor would b e consistent w i th a m ain- 
sequence-like object as proposed bv lMotch et al.l l|l987h . we 
cannot rule out the possibility that GRMus harbours an 
over-massive evolve d companion that has b een stripped of 
its outer layers fe.g. iKing fc Schenkerll2002l ). 



5 CONCLUSIONS 

We have used the Bowen emission blend and Doppler to- 
mography to trace for the first time the motion of the donor 
star in the persistent LMXB GRMus. This has allowed us 
to derive kinematical mass limits for this system. In com- 
bination with estimates for the opening angle of the accre- 
tion disc, and new measurements of the Hen A4686 emis- 
sion line (in an attempt to trace the motion of the compact 
object), we have derived tentative new mass constraints of 
1.20 < Mx/Mq < 2.64 for the compact object and 0.45 
^ M2/MQ ^ 0.85 for the companion star. 

We cannot rule out the possibility that GRMus har- 
bours an over-massive evolved companion star, or more in- 
triguingly an over-massive neutron star. Tighter constraints 
on the system parameters may therefore have important im- 
plications for our knowledge of the equation of state of nu- 
clear matter (e.g. ICook et al.l fl994) or i ndeed on the form a- 
tion scenarios for typical LMXBs (e.g. iPfahl et al.ll2003h . 

Future investigations should also concentrate on deter- 
mining the origin of the unusual second spot in the Nlll 
Doppler map (Fig. 1 11 p . This feature is hard to explain, and 
perhaps could be caused by a violation of the axioms of 
Doppler Tomography. Related to this point, the broad emis- 
sion band seen in the left hand side of the Doppler maps for 
Hell and H/3 (in this system and a number of others) is of 
uncertain origin. We believe this emission to be produced 
in an extended disc bulge, but physical modelling will be 
necessary to confirm this. 
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